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ABSTRACT 


Aircraft minimum flying speed, as determined-by actual flight test, 
is published in aircraft handbooks for pilot guidance: The test flight 
results are used to determine and confirm take-eff-and-landing speeds, 
field lengths, Teft-hand portion of the maneuvering envelopes (V-n 
diagram), etc. Determination of the absolute:minimum flying speed of 
en aircraft оп the ather- hand,, has: not been of prime.importance in 
flight test. 

Im the present analysis digital simulation’ allowed the systematic 
study of nat only the mimimum flying speed. assdefined-by Federal Aviation 
Regulations but alsa the absolute minimum flying:speed attainable in 
steady, unaccelerated flight... The: study included:such effects as 
deceleration rate, rate af change of elevator-angle;.aircraft weight 
and pitch moment of inertia.. 

It was found far an assumed light-weight-fighter aircraft that the 
absolute minimum flying: speed was. approximately 20 knots less than the 
FAR minimum flying speed.. Moreover the FAR- minimum ffying speeds 


tended to be quite sensitive to rate of change of elevator angle. 
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SYMBOLS: 


normal flight path acceleration,. ft/sec 


pitch moment of inertia,. slug-ft^ 


location of the aircraft. center of mass: 
mean aerodynamic chord 

drag coefficient, drag/qS 

lift coefficient, lift/qS 


maximum Tift coefficient 


apparent lift coefficient,. weight/qs- 


Maximum apparent lift coefficient. 


static lift coefficient as a function of angle of’attack 


rate of change of lift coefficient: for-ac change: in” 
elevator angle 


rate of change of Tift coefficient: for-azchange in: 
pttch velocity. 
pitching moment coefficient,. moment/qSc_ 
static pitching moment coefficient. as:a: function of~ 
angle of attack: 
aC. 
Eee ) 
ос: 
d+) 
pitching moment coefficient derivative for 
elevator deflection 


angle of attack damping derivative, 


oC 
=— — E 


) 


pitch damping derivative,. 


thrust coefficient, thrust/qS 
longitudinaT aeradynamic force: coefficient;.body axis 


vertical aerodynamic force coefficient; body axis: 





[L,M,N] 


[P,Q,R] 


ч YW әр e 


[u,v,w] 


product of inertia (fyzum,.slug-ft^- 


elevatar angle, deg.. 

rate af change af elevator angle, deg/sec:. 

local. truncatior error for Runge-Kutta: aigorithm 
product of inertia {kydm, slug-ft 


resultant external force: vector relativezto:aircraft 
center af mass 


resultant. external moment. vector about:the:aircraft center 
af mass 
scalar acceleration. due: to gravity,.ft/sec^ 


angular momentum vector-relative to-aireraft:center of 
mass 


scalar components: af G;: rolling, .pitching:and -yawing 
moments, ft-lb.. 


scaTar- components of'W;: rolling, pitching:and :yawing 
velocities, rad/sec.. 


rate af change af pitching velocity. 

dynamic pressure, lb/ft^ 

dynamic pressure at: minimum flying.speed/stall speed, lb/ft? 
wing area,. f 


thrust, Tb. 


scalar components ar relative to body axis system, 
ft/sec. 


rate of change of V, ft/sec? 
rate of change of Vo» as determined by FAR method 


rate of change of V_,.as determined: by altitude break 
method 2 


magnitude of resultant linear velocity of' aircraft center 
mass, ft/sec. 


resultant velocity vector of aircraft-center of mass, 
ft/sec 





вада 


[х' Y'.z'] 


гара 


тї 


Yr+1,2 


[+,0,$] 


Е | 


minimum Flying: speex? ass determined: by C,. method; .ft/sec. 
minimum Fly ing: speet ass determined: by. Cf method; .ft/sec. 


minimum Flying. speed: as: determined: by altitude:break 
method, ft/sec.. 


minimum flying speed: as: determined: by constant deceleration 
rate, ft/sec.. 


&ircraft weight,. Ib.. 
rate af change: of’ W,. scalar" component:of“V_, .ft/sec“. 


cartesian coordinate axes: notation for"body axis:system, 
argin located at. aircraft. center of mass: 


Cartesian coordinate- axes: notation for-stability axes 
system, argin located. at.aircraft.center of mass: 


components of resultant: aerodynamic: forcecacting<on 
aircraft, ТЬ... 


vartable in Runge-Kutta- algorithm obtained by intergrating 
between twa points, X, and X > with step size-h]. 


variable in Runge-Kutta: algorithm obtained:by intergrating 
between two points,. X, and.X j> with step size h, where 


angle of attack,. deg.. 
vate of change of angle of attack,.deg/sec:. 


aircraft pitch angle,. angle between horizontal reference 
and the X axis of the body axis system, deg. 


rate of change of pitch angle, deg/sec. 
air density,. sTugs/ft? 
Euler angles, rad. 


angular velacity vector of the aircraft, .rad/sec. 





L. INTRODUCTION 


Minimum flying speed and ar stalling: speed: is: defined: by Federal 
Aviation Regulations (FAR) [Ref.. 1] as: the: minimum steady speed :at * 
which the aircraft remains contrallable.. The:Navad Test:Pilot:School 
Flight Test Manual [Ref.. Z] defines: minimum flying: speed.as :the :minimum 
steady airspeed attainable in unaccelerated: flight.or"the:minimum 
usable airspeed. The actual determination of minimum flying-speed is 
always accomplished through flight. tests:. Methods: for -theoreticat 
predictions are empirical, based. upon wind-tunnel results, experience, 
ED. 

FAR specified that stalls: would. be: demonstrated:by trimming the 
aircraft at 130 percent of the: estimated minimum flying-speed -and - 
deceterating at a constant rate- until a. minimum speed.was:obtained. 

The certified stall speed would: correspond. to a: deceleration rate of 

1 knot/sec. To facilitate the test. pilot's. task,.Ref:.l and.2.defined 
aircraft characteristics and parameters: indicative:of^minimum flying 
speed. Typical indications were: longitudinal, directional or lateral 
divergence, excessive altitude loss,. loss: of control/effectiveness, etc. 
Data analysis methods and restrictions: on demonstration technique, 
other than rate of change of airspeed,. were: not: specified by Ref. 1. 

Additional flight test techniques have: been proposed and evaluated 
by the National Aeronautics and Space Administration [Refs.3 “and 4]. 
These studies, through flight tests,. evaluated three- demonstration 
techniques: FAR, l-g break and constant-rate-of-elimb. Three methods 


of data-analysis were used with the FAR demonstration technique .to 
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determine minimum flyTng speed:: constant. deceJeration;.C* and:Cr methods. 
The minimum flying speeds obtained. hy the three- techniques:were:corrected . 
to a deceleration rate of T. knat/sec im accurdance: with Ref: .1.. 

The study described herein applied. numerical methods: to: simuläte: 
flight test data in an attempt ta. define an absplute:minimum flying: 
speed. Five data-analysis methods were used.. The-speed.was:an absolute. 
minimum obtainable for steady, unaccelerated flight: when” the=FAR- 
correction to a I knot/sec. deceleration rate was: disallowed:. The: 
absolute minimum flying speeds were. compared with FAR- certified-minimum 
flying speeds, 1.e., corrected to | knat/sec.. deceleration’ rate:. 

"Flight test data” for the F-94A,. a small. single-engine: jet-fighter, . 
was obtained from a computer program of the non-4 inear’ equations :of: 
aircraft motion. FORTRAN IV language and the: Naval Postgraduate=School | 
IBM 360 digital computer were utilized.. The: computer” program provided: 
numerous options far trimming and flying the. aircraft,. including: 
variations of thrust, weight, rate af change of elevator angle_and: 
flight orientation, T.e., cTtmbing,. level or descending: ffight:. 

For this study, the aircraft was trimmed im steady,.descending: 
flight at 130 percent of the estimated minimum flying speed-and decelera- 
ted by decreasing the eTevator angle at a constant rate. Two values. 
for thrust were used to trim and stall the aircraft: T/W s.0.02_and 
T/W = 0.11. 

The effects on minimum flying speed of thrust,. aircraft:weight; 
rate of change of elevator angle and aircraft: moment: of” inertia:about 


the lateral axis were also studied. 


ТТ 





The purpose af this study. was: toc investigate-the:retationships of 
the aircraft minimum flying speeds: as: determined-by the:different data- 


analysis methods. 
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IL.. EXPERIMENTAL. METHODS: 


A. BASIC APPROACH 

The study of minimum flying speed was: restricted: to: aircraft motion 
in the XZ plane: i.e., to aircraft motion along: the: X: and-Z_axes:and. 
pitching moments about the Y axis. It was: assumed: that: no: adverse: 
lateral or directional tratts accurred while: decelerating: to: the: 
minimum flytng speed. 

Non-linear equations af aircraft motion with" three- degrees-of-- 
freedom were derived fran EuTer's equations: of" motion,. Appendix. A. . 
Extensive use was made af Etkin [Ref.. &].. The: equations: werecset:up. 
for input of non-linear aerodynamic. data,, by means: of az table: look-up. 
procedure. 

A computer program was written. im FORTRAN IV. language: for~the<Navat ` 
Postgraduate School IBM 360 digital computer,. Appendix B.. The: program 
employed a fourth-order Runge-Kutta algorithm in solving the:equations. 
of motion. Real-time was used in evaluating the: aircraft. motion’ from 


approximately 1.3 Vain to V Typical time: histories. of the 


min 
aircraft parameters were plotted in Figure Т. Figure. 2. contains summary 
plots of n as a function of deceleration rate. for different:rates 


of change of elevator angle. The thrust to weight ration (T/W) was:0.02. 
B. TEST AIRCRAFT 

The F-94A, a small single-engine jet. fighter,.was: utilized as-the. 
test aircraft. 

Reference aeradynamic data were obtained. from Blakelock [Ref:.5] in 


the form of linear stability derivatives for all terms except: those 
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Figure 1. Typical Time Histories of Aircraft Parameters. 
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depending upon angle af attack.. The variation of’the=:complete:aircraft 
с, Cp and C, was estimated from the available data: för a:center of 
gravity located at the 0.252 C (wing: meam aerodynamic:chord) station 
and an airplane in the Tanding configuration. Although the-=variations 
of C > Co and Cy with angle af attack,. a,. were estimated: for zero 
elevator deflection, the use: of the linear: stability derivatives ‘such 


as CL, and M. allowed trimming the: aircraft to any initial flight 


бе Sa 
conditions. The C, Vo ос сигуе Рог 6 — 0 was се| естед- Рог a:trim 


C of about T.46 with the stall break being very. gradual, which is 


| пах 


indicative of a progressive flow separation starting-at-the-wing trailing 
edge. 

Aircraft reference data and physical characteristics:are-listed in 
Table 1. 

The choice of atrcraft flight. traits: are arbitrary and:may seem 
restrictive to a person wha is solely concerned about:a:particulàr 
aircraft by virtue af personal experiences. However,.it-should be 
recognized that the methods and qualities. described. in" this.analysis 
are typical and may be modified to a:particular application by changing 


the input data. 


ВНЕ НТ РВОСЕВУВЕЗ 

All "test flights" were conducted with the F-94A configurated for a 
landing approach. The aircraft was trimmed for steady descending flight 
at 130 percent of the estimated minimum flying speed. Aircraft thrust, 
weight and moment of inertia about the lateral axis.were:the.only 


physical characteristics changed during the "test flights!. 
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TABLE I 


&. Physical Data 


Parameter Value 
$ 239:ft*- 
E МАЕ. 
W 12;359 1b 
D. 26;545 :slug-ft^ 


С, (Tinear region) 
а. 


зенит 


C 0643 гад! 
So | 

См (Tinear region) -01406 rad”! 
а 

с -0:88-rad | 
Se 

^ -4:268 rad”! 
Е 

e rad 
q 

С 25:24 C 


18 





b. 


TABLE L (Cont'd) 


Tabular Data (complete aircraft, о 


e e © © © © = e e e о o e e 


a 
(rad) 

.10472 
.13963 
.17453 
.20944 
.24435 
. 67925 
.31416 
.34907 
. 38397 
.41888 
.45379 
. 48869 
. 52360 
5585| 


CL 


0.552 
0.738 
2.921, 
1.093 
122 
1.346 
1.430 
1.484 
1 3843 
1.526 
1.500 
1.462 
1.418 
1..360 


Cu | 


- 0027 
-( 01 70: 
-0 0315: 
-0 0445: 
-0 0585: 
-0.0740 


-0.. 09:10: 


-0 1120: 


-4C 1320. 


-. 1520. 


-X 1720 


-0L 2075: 
-0. 2400: 


=O. 2795. 


Cp: 


0.0940: 
071150: 
0.1400. 
0.1695: 
052030: 
072490“ 
O.3120. 
0.3640. 
0.4050. 


0.4550: 


0.4930" 
0 5550. 
0 5720. 
0.6100° 





Deceleration from the trim condTtiom was: accomplished-by .deereasing 
the elevator angle, to simulate aft control stick movement:. Availàble: 
time histories, [Ref. 4[, tndicated elevator: angle: to: be: a: near linear ' 
function of real-time during the deceleration maneuver:. Therefore;.a: 
constant rate of change of elevator angle was: chosem for this:study. 
Rate of change of elevator angle was varied: between-005: and-- 1.15: 
deg/sec. for the different maneuvers.. 

The effects of thrust were. studied: by making. two: "test-flights" for 
each elevator angle schedule. Twa thrust. values: were- used: fór "the. 
"flights" which gave thrust to weight ratios: of- 0.02: and. 0.11... 

Five "flights" were made at -0.286 deg/sec..rate:of'change:of^ 
elevator angle while varying aircraft weight. to. study the-effects:of" 
aircraft weight on minimum flying speed.. The: aircraft weight:was: 
varied from 12,358 pounds ta 76,359 pounds: im increments: of~1000 pounds. 

“Flights” were made at representative rates: of change-of-elevator 
angle to study the effects of aircraft moment of inertia: about: the - 
lateral axis, on minimum flytng speed.. Values: 25 percent above<and - 
below the true aircraft moment of inertia were: evaluated. 

Data from the original "test flight" indicated the possible:exist- 
ance of an optimum elevator schedule in determining. minimum flying. 
speed. An exponential elevator schedule was. studied and.compared to 
the linear elevator schedule. The aircraft was decelerated using a 
linear schedule until the deceleration rate: fell below a:preselected 
rate. At this point the elevator schedule was: changed to an exponential 


schedule until minimum flying speed was obtained. 





D. DATA-ANALYSIS METHODS 
Five data-analysis methods were used to determine: minmum fTying: 
speed from the two main sets af "flight test" data.. These:methods were: 
1. Constant Deceleration Rate 

The linear elevator angle: schedule: resulted: in” near-constant 
deceleration rates by the test aTrcraft.. Minimum flying:speediwas: 
defined as the actual minimum speed obtained during. the:deceleration. 
For consistency, a deceleratian rate was defined as: the:slope:of'a: 
straight line drawn from Vin to pal so: This: deceleration rate:was : 
used in both the C and С methods... Plots: were. made: of" minimum flying 
speed as a function of decelerätian rate and: minimum fiying:speed.as.a- 
function of rate of change of elevator angle for the- two. values:of^ 
thrust. 

The constant deceleration rate method was: used: to: analyze-data - 
when the effects of aircraft weight and moment. of? inertià-on minimum 
flying speed were studied. 

2: C Methods 
сү was defined as the aircraft left coefficient: independent - 


of norma] flight path acceleration: 


The assumption of aircraft weight and aircraft lift being equivalent 
throughout the maneuver corresponds to considering C; as. being.an 
apparent lift coefficient rather than a true value. Minimum flying 
speed was defined by ће Сү method as the speed of the.aircraft when 


Ci was at a maximum value. Plats were made of Cr as:a: function of 
max 





V 


C" V 
max 


deceleration rate, as a function of deceleration-rate:and : UT 
as a function of rate af change of elevator angle. 
3. С Method 
Cr was defined as being the true: aircraft lift:coefficient which 


compensated for effects af normal flight: path acceleration: 


Нас. 
u f 
U 


q5: 


Minimum flying speed was: defined by the: Ci method as:the:speed .of the 
Ge 


aircraft when C, was at. à maximum value... Plots were:made :of- Бе as 


а function of deceleration rate, "C... asa function of -deceteration 
rate and D as a. function of rate of change of-elevator angle. 
4. Altitude Break Method 

The aircraft altitude. trace- had:a relatively: linear character- 
istics slope during the deceleration maneuver. With the:onset of stall 
the characteristic slope would increase significantly in:value. The 
minimum flying speed was defined as: the. speed of the.aircraft where the 
linear characteristic slope of the altitude trace could no longer be 
maintained. The characteristic slope increase was. indicative of 
increased altitude loss. 

A deceleration rate was defined as the slope of-a: straight line 


V V 


V : 
drawn from h to 1.1 is Plots were made of cn as a function 


min 
of rate of change of elevator angle. 
5. I-g Break Method 
The minimum flying speed defined by this method was the last 
speed at which the aircraft could: maintain 1-g flight:. This method 


proved unreliable in data reduction for an aircraft trimmed in steady, 
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descending flight and was thus not pursued further.. For: this: method,. 


an aircraft must be decelerated fram an initial climb condition.. 
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TEL. RESULTS AND. DISCUSSION 


A. ABSOLUTE MINIMUM FLYING SPEED 
1. Constant Deceleration Rate Method 

Тһе minimum flying speeds determined: were-plotted.:in Figure 3 
as a function of deceleration for thrust-to-weight-ratios;0.02 and 0.11. 
For a thrust-to-weight ratio af 0.02 the: absolute-minimum flying speed 
was 83.7 knots and occurred at æ deceleration rate-of-*2:60.knots/sec. 

A thrust-to-wetght ratio of 0.17 decreased: the. absolute-minimum ftying 
speed to 80.6 knots with na appreciable change. in deceleration rate. 
Figure 2 also indicated absolute minimum flying: speed: was:sensttive ‘to 
variation im optimum deceleration rate.. The: addition of~thrust reduced 
this sensitivity somewhat.. 

Figure 4 indicated absolute minimum flying: speed_was:retatively 
insensitive to rate of change af elevator.. Rates:between -0.7 and -l.] 
deg/sec. were able to give a« close: approximation-of absolute minimum 
flying speed. 


2. C Method 


Maximum values of C were plotted in Figure.5:as a function of 
deceleration rate. lt. was plotted as a function of-deceleration rate 
in Figure 6 and as a function of rate of change of elevator angle in 
Figure 7. A С of 1.815 was the maximum attainable in steady, 
unaccelerated flight. The absolute minimum flying-speed-was 83.5 knots. 
The thrust-to-weight ratio of 0.11 decreased: the. absolute minimum flying 
V 


C'. was 


speed three knots with minimal effects: on deceleration-rate. . An 


fairly sensitive to deceleration rate and. insensitive to. rate of change 


of elevator angle. 
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у, kts/sec. 


2. Т/М < 0.02 
Er T/W — 0.1T 





- 70 80 90 100 110° 120 
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Figure 3. Minimum Flying Speed as a Function of 
Deceleration Rate. 
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Figure 4. 
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Minimum FTying Speed as a Function of-Rate of 
Change of Elevator Angle. 
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ү, kts/sec. 


Figure 5. 





C! 


тах 


Apparent Maximum Lift Coefficient as a Function 
of Deceleration Rate. 
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у, kts/sec. 





p 80 90 100 110 


Figure 6. Minimum Flying Speed as a Function of 
Deceleration Rate. 
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JE deg/sec. 





80 90 100 110 


Figure 7. Minimum Flying Speed as a Function of Rate 
of Change of Elevator Angle. 
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2. С Method 


Absolute minimum flying speed was 88 knots far & thrust=to=- 
weight ratio of 0.02, Figures 8 and 9. Increasing thrust. decreased: 
this speed 2.5 knots. "Ga Was sensitive ta bath deceleratiom rate: and: 
rate of change of elevator angie. Figure 10 indicated the: absolute: 
max imum CL value was insensitive to deceleration rate and: increased: 
thrust. 

4. Altitude Break Method 
n was plotted as a function. of deceleration rate: in Figure: 
11 and as a function of rate of change af elevator angle: im Figures 122. 
The absolute minimum flying speed for a thrust-to-weight: ratio: of- 0:02. 
was 83.5 knots at a deceleration rate af 2.6 knats/sec.. Figure: 11 
indicated that two values for ob. existed for one value: of deceleration 
rate between 2.6 and-3.76 kts/sec. Figure TZ indicated: 'h ip wass 
insensitive to rate of change of elevator angle. 

The altitude break method. was dependent on graphical interpreta:- 
tion for data. This accounts for the smaTT amount of data: scatter im 
Figures 1] and 12. 

5. Comparison of Methods 

The general relationships of absolute minimum flying speed as. 
determined by the various methods are indicated in Figure 13. 

The constant deceleration, CL and altitude break methods defined 
the absolute minimum flying speed to be 83 knots, Figure 13. The €; 
method absolute minimum flying speed was 88 knots. Increased. thrust: 


decreased the absolute minimum flying speed defined: by the four methods. 


about three knots. 
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Figure 8. Minimum Flying Speed as a Function of Deceleration 
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Figure 9. Minimum Flying Speed as a Function of Rate of 
Change of Elevator Angle. 





ү, kts/sec. 





True Maximum Lift Coefficient as a Function 


Figure 10. 
of Deceleration Rate. 


33 





V, , Kts/sec. 


Figure 11. 


Ф. Т/К: < 0.02 


"En. T/W = 0.11 





90 100 TO 


Minimum Flying Speed as a Function of Deceleration 
Rate. 


34 





бе’ deg/sec. 





| | | 2. ТИМ 20202: 
LEE. T7W. 20:11. 
1.0 
0.8 | 
0.6 | 
| 
0.4 | 
0.2 | 
а 
80 90 100 110 
V s KES. 
m 


Figure 12. Minimum Flying Speed as a Function of Rate 
of Change of Elevator Angle. 
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The deceleration rate was approximately.-2. 6: knots7sec:at:the: 
absolute minimum flying speed far all four methods:.. Significant:changes: 
in deceleration rate due to increased thrust was: observed-only in-the: 
altitude break method. 

Minimum flying speed determined by the C, method: occurred -at - 
the smallest rate of change of elevator angle,. -0.4 degs/sec:. Minimum 


flying speed determined by the constant deceleration,. Cy  and:altitude: 


ins 
break methods occurred at -0.9 deg/sec.. rate of change: of elevator angte:z 
Absolute minimum flying speed deftned by the C, method. occurred at -the: 
smallest elevator angle, approximately -9 degrees,. whereas: approximately 
-13 degrees elevator angle was required for the other methods:. 

The angle of attack at the absolute minimum flying: speed -was - 
approximately 28 degrees for the constant deceleration,.Cy and.altitude : 
break methods, and constderably Tess at 21 degrees: for the:C, method. . 
Thrust effects on angle of attack at nh absolute: minimum flying.speed: 
were insignificant (less than 0.5 degrees) for all methods... 

Elevator angles, lift coefficients and. angle: of" attacks:at- 
minimum flying speed/stalT speed were less than,. in some:cases;.the: 
maximum static values defined for the aircraft.  Possible:explanations. 
are: 

1. Rate of change of elevator angle was. programmed: to: become 
zero two seconds after n was obtained. For small .values:of-'etevator 
rate, im occurred much sooner than the maximum obtained values:of °C, , 

a and бе. Therefore, elevator angle was held constant until the:end of 
the maneuver, restricting the aircraft from obtaining the maximum values 


of С, a and $e" 





2. The variations of drag coefficient with angle of attack: used: 
in the table look up were representative values and not: necessarily true: 
values for the F-94A. Reduced drag coefficients would: thereforezeffect: 
the occurrence of minimum flying speed. 

3. In defense of the Tow values of C >, e and: 8, ,- 1t:must:be: 
remembered that this is a dynamic study af stall traits.. Actual flight. 
test data, [Ref. 4], indicates angle of attack is: less: than" static:stall' 


angle of attack. 


B. FAR MINIMUM FLYING SPEED 

The minimum flying speed defined by Federal Aviation Regulations: 
corresponds to the minimum speed obtained at @ one knot/sec:.decelera-- 
tion rate. The FAR minimum flying speed determined by the. constant: 
deceleration, CL and С methods was approximately 105 knots. The: 
altitude break method's FAR minimum flying speed. was: 109: knots,. Increased: 
thrust decreased the altitude break minimum speed two. knots: and decreased: 
the minimum speed defined by the remaining methods: less: than one: knot. 

The rate of change of elevator was - 0.08 deg/sec..at:minimum flying. 
speed determined by the constant deceleration rate, C, and. Cy methods. 
The altitude break method rate of change of elevator angle.was -0.05 
deg/sec. 

Angle of attack was approximately T5 degrees and elevator angle -4.8 


degrees for all methods. 


C. COMPARISON OF MINIMUM FLYING SPEEDS 
The absolute minimum flying speeds were categorically 20 knots: less 
than those defined by FAR one knot/sec. deceleration rate,.Figure.M. 


Thrust had more effect on the absoTute minimum flying speed whereas rate 


i 
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Figure 14. Comparison of Aircraft Parameters at FAR Minimum 
Flying Speed. 
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of change of elevator angle had more effect an the FAR minimum flying: 
speeds. Angle of attack and elevator angle were significantly greater 


with the absolute minimum flying speeds. 


SENSITIVITY STUDIES 
1. Aircraft Weight 
For a representative rate of change of elevator angle,. Figure: 
15 indicates minimum flying speed was a linear function of aircraft,. 
weight. 


From the definition: 





minimum flying speed is seen to be a function af the one-half power-of- 
weight. Actual flight test data indicates minimum flying speed is: a: 
linear function of EU weight [Ref. 4]. 
2. Aircraft Pitch Moment of Inertia 
When the aircraft pitch moment of inertia was increased or: 
decreased 25 percent, there were no significant effects on minimum 
flying speeds or other aircraft parameters. This conclusion was. not 
based on a systematic study but rather a brief Took at three data 
runs. The trends were present but a complete study is warranted. 
3. Exponential Elevator Rate 
Comparative plots of linear and Tinear plus exponential elevator 
rates were made in Figure 16 for initial Tinear rates of -0.06 and 
-0.115 deg/sec. The exponential rate was e* where x was increased 0.04 
each second. V ., obtained by the linear plus exponential elevator rate: 


min 


was significantly less than that obtained by the Tinear elevator rate. 
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These results, though only. a sampling: of- data:,. tend:to support the, 
existance of an aptimum elevator schedule in determining minimum flying 


speed. 
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IV. CONCLUSIONS AND. RECOMMENDATIONS: 


Digital simulatian has great potential im the: evaduation-of minimum 
flying speed. The options available are: practically. unlimited: 

The determination of absolute minimum flying: speed:and:aircraft: 
behavior with digital simulation will enable the. engineer to-better 
predict the aircraft characteristics encountered: during: test-flight: 
The test pilot's job wiIT be simplified. somewhat: if^ hec has: knowledge :of 
general aircraft characteristics prior ta flight.. 

The operational pilot's performance should be: enhanced-:by having. 
an insight into his aircraft's performance. beyond: the: certified. stall. 
speed. 

The computer program coded for this study may be: used. för many 
fields of study as indicated by Ref.. 7,, where it. was: used to:make:a: 


stability analysis in the vicinity af minimum flying: speed: . 
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APPENDIX A 
NON-LINEAR. EQUATIONS OF AIRCRAFT: MOTION 


I. Euler's Equations af Aircraft Motion. 

The aircraft was considered. as e rigid: body with a: Cartesian 
coordinate system fixed at the center of mass,. Figurecl7.. Thecframe-of 
reference Cxyz, a body axis system,. moved with the: aircraft <and:Cx was 
fixed to the longitudinal. axis af the aircraft.. The:resultant external 
forces and the moments af these: forces: about: the: aircraft:center of mass, 


when referred to the frame of reference Cxyz;,.where: the:vector equations 


of motion: РЧ 
2o 3]y _ 
Реп Зпах 
БАЗА oxh 


The scalar components of these eguations. are: Euler's:equations:of- 


aircraft motion. 


II. Aircraft Orientation and Axes. 

Reference frame Cxyz could not be used. to describe the position and 
orientation of the aircraft as it was fixed. to the center of mass :and 
moved with the aircraft. Ап earth fixed: reference frame, assumming 
the earth's rotation was negligible,. was: chosen to describe the aircraft 
motion. 

Euler's equations of motion are valid. for any orthogonal reference 
frame fixed to the aircraft center af mass:.. To: simplify the-equations 


of motion and expressions for aerodynamic: fórces,.the:reféerence frame, 
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Cx'y'z', was chosen as stability axes, Y.e., Cx' pointed: in^ thecsdirection- 


of motion of the aircraft center af mass, Figure 17.. 


III. Assumptions 

In addition to the assumpttans previously ‘made: defining: aircraft- 
orientation and axes it was assumed. that aircraft. motion’ was: restricted: 
to the Cxy, longitudinal plane. This assumptiom equated: the: following: 
quantities to zero: L, N, P, R, &, w, V and Y... X Cxy. was: also: assumed: to: 
be a plane of symmetry which equated the products: of. inertia;.D.and.F;,. 


to zero. 


IV. Non-Dimensitonaltzation 

The external aerodynamic forces and moments: were: expressed_as: 
functions of dimensionless force and moment. coefficients. .dynamic<pressure, 
aircraft planform wing area and a reference: length.. Thesecequations;. 


in scalar form, were: 


5; 


>< 
И 

N 
n 
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> 
AN ARN 
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Г ое ее 
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When these values were substituted inta the scalar resultant external 


force and moment equations the following equations were obtained: 


| gest, 

[IQ a sin 6 — QW 
2m 

rc, 

W = + g œs o + QU 
2m | 
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Figure 17. Aircraft Axes. 
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V. Aerodynamic Force and Moment Coefficients: 


From Figure 18, the force coefficients: were: 

on ee іП — C RO 
т C G sin a. Cy COS: a: 
- er cos a — Cy sin a 


The pitching moment and [Tft coefficients: consisted of-linear portions. 


N 


C 


" 


The non-linearity was present in Q (a) and' C (a) ,. which were:non-linear 
curves of С апа Cn versus . Complete expressions: fór-the:total' 


pitching moment and [ift coefficients were: 


єз 
| 


е. е | 
j \ ( k - ^ = + 
- C(ye&£lc 6+ + | 
m Cfo) дү, | | Cir Ё | | С: | бе: 
q аз Ji Ša: 


*" = Вазов езс в. 
L Gi ke dE: Lse- = 


VI. Final Equations 

Angle of attack appeared in the pitching: moment equation as a 
derivative. To obtain a set of simultaneous differential equations that 
could be solved by one of the methods available,.another equation was 
added defining the rate of change af angle of attack. The three-degrees- 
of-freedom, non-linear differential equations: of aircraft motion programmed 


for the Naval Postgraduate School IBM 360 digital computer were: 
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Figure 18. Aerodynamic Force and Moment Coefficients. 
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APPENDIX B: 
DEVELOPMENT OF- COMPUTER- PROGRAM 


I. COMPUTER FACILITY 

The ariginal intent. was. to program the:equations.:of motion on the 
Aeronautical Engineering [leepartment: EAI.580.analog computer. Problems 
with insufficient equipment to. program thesequations :of motion on the 
enalog computer necessitated the:use-of' digital simulation. 

The Naval Postgraduate: School's: IBM 360-digitat computer offered 
many advantages and was: chosen as: the: research tool... Several major 
advantages were:. the large: amount: of corezavailable,; .dimensioned 
quantities could. be used,. plot. subroutines:were available «for time 
history studies and. numerous methods werezsavailable.for solving 


simultaneaus non-linear: differential equations on a :digital computer. 


II. COMPUTER PROGRAM 

ОҒ the many one step procedures: that:were:avatlable for solving 
systems of simultaneous: first-order- ordinary. differential equations, 
a fourth-arder Runge-Kutta algorithm was:chosen, [Ref: 8]. This method 
offered simplicity, accuracy and: short computer time. The Runge-Kutta 
procedures were programmed: into subroutine RUNKUT.  RUNKUT called 
functicn subroutines Fl, F2, F3, F4 and F5 to calculate the derivative 
terms at time T(I), T(I)-HI/2 and. T(I)4HI. The:derivative terms were 
multiplied by the appropriate. weight. factors and.added to the variables 
calculated at time TI) which gave. values.for the.variables at time 


T(I+T). 
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Function subrautine FT, c&lleıt subroutine: SPLIN to-calculäte:the: 
non-linearities în G., C, and C... SPLIN. wass obtained: from the<sctentific 
subroutine package af the Naval Postgraduate: School IBM 360-digitat 
computer and modified for this program.. Non-linear tables:of C.,.C. and 
Cp as a function of angle af attack were read: into: SPLIN.. When "supplied 
values of angle af attack, SPLIN. used. @ cubit: spline: function: ta provide 
interpolated values af Q» Cn and C. 

Initial aircraft trim conditions U,. W,.0,.6,.a.and-6, were supplied 
to RUNKUT at time zero far initiation of the:step:procedure.. These: 
values were calculated in Subroutine. TRIM from static-fórcecand moment 
equations. 

Control of the aircraft was accomplished: through changes:in elevator 
angle and thrust which simuTated control stick and.throttle:movements. 
The options available were selected through the: integers:on” 


the first data card. 


Ill. TRUNCATION ERROR AND STEP-SIZE ANALYSIS. 
Error analysis by rigorous mathematical derivations:are:virtuatly 
impossible to implement for higher-order Runge-Kutta algorithms -for 


systems of differential equations. The truncation error: 


—À 


6 


"t 7 Sox ry 


for a first-order ordinary differential equation was:applied to the 
system of differential equations as a step-size control mechanism, [Ref. 8]. 


Step-size and truncation errors are listed in Table. II. 





Table II 


Local TruncatTon Errors 


Step-Size Error 


0.2 -Z.92 x W^ 
0.1 2.97 x 107° 
0.05 1.55 xx 107 9 
0.025 ETE ES 
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